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We study the stability of AdS black hole holes rotating in a single two plane for tensor-type 
gravitational perturbations in D > 6 space-time dimensions. First, by an analytic method, we show 
that there exists no unstable mode when the magnitude a of the angular momentum is smaller 
than r\/R where is the horizon radius, and R is the AdS curvature radius. Then, by numerical 
calculations of quasinormal modes, using the separability of the relevant perturbation equations, we 
show that an instability occurs for rapidly rotating black holes with a > r\/ R, although the growth 
rate is tiny (of order 10~^^ of the inverse horizon radius). We give numerical evidences indicating 
that this instability is caused by superradiance. 



I. INTRODUCTION 

Recent years, the perturbative study of higher-dimensional black holes became an active research area P, 0| along 
with the developments in string theory requiring extra dimensions and various higher-dimensional world models 
iii. 

In such a world model, the black hole stability is a crucial problem because no simple uniqueness theorem holds in 
spacetimes with dimensions higher than four. Instead, it has turned out that a great variety of higher-dimensional 
objects can exist in higher dimensions: Myers-Perry (MP) black holes, Kaluza-Klein black holes, different types of 
black strings, branes, rings and saturns [lO|. In order to know which of these solutions can be realized in nature, we 
need to test their stability against small metric perturbations. The black hole stability is also conjectured to provide 
useful information on the phase structure of some conformal field theories through the AdS/CFT correspondence 

sag. rnr:.r-^ 

The systematic investigation of the black hole stability in higher dimensions started from the papers [ll|, Il3 ' 

where the master wave equations for gravitational perturbations of Tangherlini-like static black holes allowing for the 
cosmological constant A were derived. Then, with the help of them, the stability of higher-dimensional Schwarzschild 
black holes was proven analytically [3]. The analytical tools used there, however, did not work well in the presence of 
the non- vanishing A-term or electromagnetic charge. Therefore, the stability problem in these cases have been studied 
numerically. For example, the stability of D-dimensional Schwarzschild-de Sitter black holes was shown numerically 
in Ti"], and that of Reissner-Nordstrom (RN) black holes in [Ts']. Further, recently, a gravitational instability has been 
found for RNdS black holes in Z? > 7 space-time dimensions, when both charge and A-term are large enough ^15 ]. 

When we extend the analysis to the system with negative A and non-vanishing charge, a subtlety arises. It is 
because the Reissner-Nordstrom-anti-de Sitter (RNAdS) metric can be a solution in the D-dimensional Einstein- 
Maxwell theory and in the N — ?> gauged supergravity theory at the same time. The latter theory includes a dilaton. 
Thus, the dynamics of black hole perturbations evidently depends on the theory in which we consider the RNAdS 
metric. For example, the RNAdS black holes are dynamically unstable in the M — % gauged supergravity with I? = 4 
when the AdS black hole is large and highly charged [7], while the RNAdS black holes in the pure Einstein-Maxwell 
theory are dynamically stable [8']. 

Motivated by the low-energy limit of string theory, D-dimensional black holes in the Gauss-Bonnet theory were also 
studied and proved to be unstable for large values of the Gauss-Bonnet coupling and only for Z? = 5, 6 space-time 
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dimensions [la. llTl . The stability and quasinormal modes of Kaluza-Klein black holes with squashed horizon was 
considered in jl8l|, while the evolution of instability in time domain for black strings was analyzed in 19]. 

Thus, the stability of higher-dimensional static black holes has been studied quite well up to now (though not 
completely). Unlike the static case, the stability of rotating black holes in higher dimensions has not been proven so 
far, because the perturbation equations for them cannot be reduced to a set of decoupled wave- like equations and the 
separation of variables is a hard problem as well. Nevertheless, many studies suggest that rotating black holes tend 
to be more unstable compared with static ones. For example, when we consider a bosonic field around a rotating 
black hole, there appear modes with frequencies Rew < mflh that extract a part of the rotational energy of the black 
hole at scattering. For such a mode, a reflected wave has a larger amplitude than an incident one. This effect of 
amplification of the incident wave was called the superradiance and calculated for the first time for Kerr black holes 
by Starobinsky [2l[. Because the effective potential for a massive field has a local minimum, the wave amplified by 
superscattering can be reflected back onto the black hole again by the wall around the minimum. The repetition of 
this superradiance and reflection can lead to the growth of the wave amplitude without bound for a massive field 
around a Kerr black hole [IHl , although its growth rate is in most cases too small to have significant effects in realistic 
situations. 

This type of instability is also expected to occur for rotating AdS black holes because the asymptotic behavior of 
the AdS spacetime produces a diverging effective potential barrier at infinity even for massless fields [2^. In fact, 
rotating AdS black holes with all equal angular momentum components in spacetimes of odd dimensions with D > 7 
were shown to be unstable for a special class of tensor- type gravitational perturbations in (20| . 

The main purpose of the present paper is to investigate the possibility of this type of instability for another special 
class of black holes that are much more physically motivated: the asymptotically AdS simply rotating black hole, 
which is a black hole with only one non- vanishing angular momentum component. Indeed, if black holes are created 
by particle collisions, the colliding particles, being located on our four-dimensional brane, will create a black hole 
rotating in a two plane contained in the brane dimensions. One specific feature of such a simply rotating black hole 
is that the metric has a sphere component and as a consequence, its gravitational perturbations can be classified into 
scalar, vector and tensor types and expanded by tensor harmonics as in the case of static black holes. In contrast 
to that case, however, the perturbed Einstein equations cannot be reduced to ordinary differential equations (ODEs) 
by this harmonic expansion even after the Fourier transformation with respect to time and the rotational angular 
coordinate. Nevertheless, the problem can be reduced to a single two-dimensional partial differential equation (PDF) 
for tensor-type gravitational perturbations, which turns out to be identical to the massless free scalar field equation 
on the black hole background and can be reduced to ODEs by the separation of variables. 

In the present paper, we use these features to study the stability of simply rotating higher-dimensional AdS black 
holes for tensor-type gravitational perturbations. Because the tensor-type perturbations on S" exist only for n > 3, 
the analysis applies only to the spacetime dimension D = n + A > 7. We investigate the stability by two methods. One 
is an analytic method motivated by the general arguments of Hawking and Reall [1^ . By this method, we prove that 
the simply rotating black hole is stable against tensor-type gravitational perturbations in the slow rotation regime 
a < rf^/ R, where rh is the horizon radius, and R is the AdS curvature radius. This method, however, cannot be used 
to find an instability. Hence, we perform a detailed numerical analysis of the quasi-normal frequencies for tensor- 
type perturbations. We will show that in addition to the decaying part of the spectrum corresponding to damped 
quasinormal modes, there are unstable modes that can be interpreted as being produced by superradiance. Indeed, 
we numerically confirm with very good accuracy that superradiant modes are unstable and that non-superradiant 
modes are stable in all the cases we calculated. 

The paper is organized as follows. Section II introduces the basic features of the background metric for the simply 
rotating black hole. In Sec. Ill, we give the tensor-type perturbation equations for this black hole and show that 
it can be reduced to two weakly coupled wave-like equations. In Sec IV, the stability in the slow rotating case is 
proved analytically with the help of an energy integral. Section V explains the numerical methods to determine the 
quasi-normal frequencies by solving the weakly coupled mode equations for the angular part and the radial part. 
Section VI discusses the numerical results obtained for damped quasi-normal and unstable superradiant modes. In 
Sec. VII, we summarize the results and comment on open questions. 

II. BACKGROUND SOLUTION 

In this section, we summarize the basic features of the simply rotating asymptotically AdS black hole solution in 
higher dimensions that is used as a background for perturbations. 
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A. Metric 



The general higher-dimensional rotating black hole solution with non-vanishing cosmological constant found by 
Gibbons, Lu, Page and Pope [23 | has A^-independent angular- momentum parameters ai, • • • ,ajv, in addition to the 
mass parameter M, where N = [{D — l)/2] for Z?-dimensional spacetimes. In the asymptotically flat case, i.e., for the 
Myers-Perry solution, these parameters correspond to the eigenvalues of the antisymmetric matrix Jab representing 
the total angular momentum, and each of them to the rotation in one of the N mutually orthogonal 2-planes p3| . 

Although the expression for the metric for generic momentum is rather complicated, it becomes quite simple for 
the simply rotating case in which the black hole rotates only in a single two plane [l^l : 



ai — a, a2 = ■ ■ ■ = cin = 0. 



(1) 



The metric for the higher-dimensional AdS-Kerr black hole in this case can be written in the Boyer-Lindquist coor- 
dinates as 



X 


■ 2M 
















X 



-X - Xo{l-\r'^)p^ 



, 2 AaMX sin^ 6 , , 



xIp^ 



Xa{r'+a')p' + 



2a^M 



COS^ QdVl} 



where dfi^ is the metric of the n-dimensional unit sphere, and 



A:= (r2+a^)(l-Ar2) 



2A/ 



2 2 

p := r 



= 1 + \a\ 



(2) 



(3) 
(4) 



cos^ 6*, X = 1 + Aa^ cos^ 

This metric can be written as the direct sum of the metrics of a 4-dimensional base space and spherical fibres as 

ds" = 5afc(2/)dy"d2/' + S{yfdVL\,_^, (5) 

where a,b = t, 0, r, 9 and 

5'(y) = rcos(6'). (6) 

The determinant of the base metric is given by 

det(g„b)i/2 = |l-FAa2|-V^sin6l. (7) 



From the above or the expression ([2]), the metric is singular for Xq = 1 + Aa^ = 0. However, it is not obvious from 
these expressions whether Xq < is allowed or not. In order to see that this parameter range is not physical, let us 
rewrite the metric in the form 



,9 ^ / , a . n 

ds^ = 2[dt- sm^ 



X sin^ 



adt — 



C 



'--dr'^ + —d9'^+ cos^ Odn^ 



X 



D-4- 



(8) 



If Xq < 0, the function X changes the sign when 9 varies from to tt/2. However, from this expression, it follows 

that the metric has the wrong signature [ h • • • +] for X < 0. Hence, the metric ([2]) is regular for all values of 9 

only when Xq > 0, i.e., < — 1/A. 



B. Horizon 



From 

9t<i, - 9tt94>4> = sin^ 9, 

the location of the Killing horizon can be determined by the equation 

2M 



(9) 



(10) 
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Since the spacetime is singular at r = 0, in the five dimensional case D = 5, the black hole has a regular horizon 
only when 

< 2M. (11) 

Hence, there exists an upper bound on the angular momentum parameter a, as in the case of four dimensions. In 
contrast, for D > 5, A = has always a single positive root if A < and M > 0. 
The horizon angular velocity flh is determined by the condition 

k = dt + ntd^, k-k = <^ nlg^^ + 2nhg^t + gtt = 0, (12) 



as 



_ a{l-Xrl) _ 2aM 



III. TENSOR PERTURBATIONS 

For spacetime dimensions greater than or equal to five, from the metric structure we can decompose gravitational 
perturbations of the simply rotating Kerr-A black hole into reduced tensorial types according to the tensorial behavior 
with respect to S^~'^ [2l |25||. In particular, for D >7, tensor- type perturbations appear. 

A. Master equation 

In general, tensor-type metric perturbations for the background ([5]) can be expanded in terms of a basis of transverse 
and traceless harmonic tensors on the unit sphere S^^'^, as Q 

Sgab^O, <5g„ = 0, - 25(y)2 ^ 4''")(y)T5"'(z). (14) 

Here, i and j refer to the coordinates z — (z*) of the sphere S^~'^, I labels the eigenvalue of the harmonic tensor 



^'^^ ^^'^ Laplacc-Bcltrami operator A on S as 



(I 



A + l{l + D-5)-2 T^y ^ =0, (15) 

and a is the label to distinguish harmonic tensors with the same eigenvalue. 

After this harmonic expansion, the Einstein equations can be reduced to decoupled second-order hyperbolic equa- 
tions for each amplitude h!^'°'\ Omitting the index {I, a) for simphcity, this master equation reads [^|25| 

- DHt ~ ^g^'daSd,HT + ^-^^^^^Ht = 0, (16) 

where □ denotes the d'Alembertian for the metric gab{y)- Note that this equation is completely identical to the 
equation obtained from the free massless field equation for a scalar field on the same background by the harmonic 
expansion with respect to the harmonic functions § satisfying 

A + l{l + D -5)\s = 0. (17) 

After Fourier expanding Ht with respect to t and as 

Ht = H{r, 9) exp(-iwi + im(j)), (18) 
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where m is an arbitrary integer, we obtain 



-D + 2- ^ + A {Dr^ + {D- 2)a^} 



rD-3 



rdrH 



'{l + Xa^ cos^ e)d^H + 



-- il + ^a^){l-Xr')- 



{D - 4) sin^ e ~ cos2 
sin 9 cos 9 
2\M \ l + \a^ 



{D -2- [D -I) cos^ 6*} cot ( 



A 

+Z(/ + -D-5) ( ^ 



r-D-5 



sm 



2 • 2 , 

a sm I 



(r2 " 


cos^ 9 ) 








B. 



1 + Aa^ cos^ I 
H = 0. 



B. Separation of variables 



(19) 



The equation for H can be separated in the foUowing way 

H = P{r)Q{9) 

to the angular part 

' (D - 4) sin^ 9 - cos^ a 



(1 + cos^ 9)Q" 
1 + Aa2 



Aa^ {£> - 2 - (£> - 1) cos^ 9] cot 6* Q' 



9 • 2 

a" sm I 



and to the radial part 
2M 



sin^9 



1 + Xa? cos2 6' 



cos^ 



AP" 



r 



D-3 



D - 2+^-^^^^ - X{Dr^ + (D ~2)a^} 



rP' 

,2 



^ , (1 + Aa )(1-Ar )--^ I --^^amtu+ 



(20) 



(21) 



P = 0, (22) 



where fi is the separation constant. 

Thus, the stability analysis is reduced to the eigenvalue problem of ODEs as in the case of static black holes. In the 
rotating case, however, there is a subtlety. In contrast to the static case, oj^ also appears in the ODE for the angular 
part. Hence, these two equations are coupled as the eigenvalue equation. 



IV. ANALYTIC STUDY 



In this section, we prove the stability of black holes with A < — against tensor perturbations. The basic idea is 
to apply Hawking and Reall's general argument 26] for the stability of a scalar field in a slowly rotating Kerr-AdS 
background to the present system that is equivalent to such a system. Although it is sufficient for the stability proof 
if we show that there exists an everywhere causal Killing vector outside the horizon (27'|, in the present paper, we 
explicitly construct the energy integral corresponding to that Killing vector, because its structure gives us information 
on instability in the rapidly rotating case studied numerically in the next section. 

In the stability argument in the present section, we do not utilize the separability of the differential equation for 
H(r,9). Instead, we directly work with the function $ of r and x = 008(26*) defined by 



$ obeys the following equation: 



H 



-£1/2+2/2 



(^2 + a2)-i/2(i + a;)-(^-5)/4a>(r,a;). 



(23) 



.F 



{lo — mf2) — + m(a; — mfi) 



{2(1 - x^){2 + Xa^ + Xa^x)d^<i>} 
4a^{l - x)M 



(2 + Aa2 + Aa2x)(r2 + a'^^r 



2^2„_D-5 



+ m^Ui{r,x) + Uo{r) 



$ = 0, 



(24) 



where 



_ 2aM 

_ (1 + Aa^)(r^ + a^)(2r^ + + a^x) + 

(£> - 2)^4 + 2(13^ - eJ) + 4)aV^ + (£> - 4)^a^ 
° 2(r2 +a2)2r^-3 

3 (iJ-4)(i? + 2)(l + a.) + 7(l-a:) , 2 ^(^-2) , 2 

+ 4 8 

a2(l + Aa2) /, \f, D \ {21 - b + Df 



Ui := 



(l + Aa2)(2r2+a2 + a2x) 2a2M 



(1 - a;)(r2 + a2) (r2 + a2) 

(2 + Aa2 + Aa2x)(r2 + a2)4^2D-io ' >^ °) 



A. Effective potential 

In the static background, the eigen-frequency uj for an unstable mode is purely imaginary, that is, < 0, because 
the corresponding differential operator is self-adjoint. In contrast, the superradiance arguments [27] and the stability 
analysis for the special Kerr-AdS black holes 20] suggest that an instability mode appears at the complex frequency 
whose real part satisfies the superradiance condition lo — mflh = 0. This is natural because defined by 



Lo — mflh 



(29) 



is the frequency with respect to the advanced time on the future horizon that is regular with respect to the null 
geodesic generator of the horizon [131 ■ 

On the basis of this observation, we take UJ^, as an effective frequency. In terms of this, the equation for <I> can be 
rewritten as 



(r^ + a^)dr 



A 



{2(1 - x^){2 + Xa^ + Xa^x)d:r'^] 



-uj;— - 2mauj^G + m^'Ui + C/q $ = 0, 



(30) 



where 



a 



(l + Aa2)(2r2 + a2 + Q^x) 
(1 - Ar2)(r2 + o2)(2 + Aa2 + Aa2a;) 

1 ^ Xr^ + 5{l + Xa^){r^ + a^) + 

{l + Xa^f{2r-'+a'+a'xf 
(1 - x)(r2 + a2)2(2 + Aa2 + Xa'^x) 



2da^M 



1-x 



2r2 + a2 + a'^x 



1 - 



^2(r2 + a^){l - x) {2,5(1 - Ar^) + ^2(^.2 + a^){i + Xa^)} 



(1 - Ar2)(2r2 + a'^ + a?x) 



2a'^M5{l - x) 



^D-5(i _ Ar2)(2r2 + + a^x) 



{2 + 5(r^ + a') 



, 2(1 + Aa^)(r^ + a^) -a^jl^ x){2 + Xa^ - Ar^) ■ 
(1 - Ar2)(2r2 + a2 + a2.T) 



(31) 
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In these expressions, ^ is a function of r defined by 



2 „2 



(5 = (1 + Aa^)^ 



(33) 



which is always positive outside the horizon. 



B. Behavior of the potential 

G is always positive, and f/i is also positive at the horizon r = rh- 

(l + Aa2)(2r| + a2 + a2a;)2 



Uiirh) 



(1 - x)(a2 + rl){2 + Aa2 + Xq-^x) 
However, for A < 0, f/i at r — > oo can become negative as 



> 0. 



4(1 + Aa2)2 
2 + Aa^ + Aa^a 



1 + a; 



+ Xr 



2{l-x) 2{-X){rl + a^r 



(34) 



(35) 



Now, wc show that C/i is non- negative everywhere outside the horizon if this asymptotic value is non-negative. For 
that purpose, we introduce the parameter b by 



|A|' 



which is in the range 

< 6< 1. 

In terms of b, U\ can be expressed in the units r/j = 1 as 

{l-bf 



(36) 



(37) 



(1 - x)|A|(6 + |A|)2(2 - b - bx)r^-^A 



2 r \ 2/r,D-2 „D-4 



b^X\2' 



l) + {b+ l)|A|r(r^-3 - 1) + 6(r^-4 - r)} (1 + xf 



+2b\X\ {2AV(r^-i - 1) + (6 + l)|A|(r^-3 + r'^-^ - 2r) + 2b{r'^-^ - r)} (1 + x) 
+4A2(&|A| +2b+ |A|)r(r^-3 - 1) + ^X^{X^ - b)r{r^-^ - 1) . 



(38) 



All terms except for the last one in the above expression are non-negative. Further, for a; = —1 adn r ~ oo, the last 
term dominates. Therefore, the necessary and sufficient condition for f/i > is expressed in the general units as 



X^rl - 6 > <^ < -Art 



(39) 



C. Energy integral 



Suppose that ^{t,r,x) is a solution to the equation 

A 



-(r^ + a'')dr I ^ dr^ I - {2(1 - a;2)(2 + Aa^ + Aa^a;)^^*} 



+ -rd^^ - 2imaGdt^ + {rn^Ui + Uq)^ = 0. 



Then, for the energy integral defined by 

J^($) := / fir / dx 
F 

+ 



(r2 + a2)A 



+ a 

|a,$|2 + !!!!^±J^ 11.12 



(40) 



(41) 



8 



pure adS limit 




small bh-* 



large bh 



FIG. 1: The stable region in the parameter plane for the simply rotating MPadS black hole 



we have 



dr / dx 



-drdt^*dr^ + 2^ \ . ^ '-d.j:dt^*d^^ 



dx 



F 



, (r2 +a2)A 



CC 



dr \ dx- 



A 



-d.^ |2(1 - x^){2 + Aa^ + \a^x)d^^^ + ^8"^^ + {m^tji + C/o)«> 



CC 



, , , 2imaG 

dr I dx— -\dt<^\^^cc 

_i r^ + 



= 0, 



(42) 



where we have assumed that $ falls off sufficiently rapidly at r ~ cx) so that is finite. 

Hence, Jff{^) is constant in time. From this and Uq > 0, it follows that if f/i > holds, (00]) does not have a 
normalisable solution of the form $ — $(r, a;)e~'"** (Im > 0). In particular, the black hole is stable against tensor 
perturbations if A < —a^/rf^, or equivalently a < r^/R, where A — (see Fig. [1]). Taking into account the 
regularity condition a < R, this condition can be rewritten as Rflh < 1 because we have 



iR-a){aR-rl 

HI i/i - i — ^ — — — 

R(rf^ + a^) 



(43) 



from (Uni). 

Finally, note that the arguments in this section also imply that the stability of modes do not depend on the sign of 
m and that instability can occur only for a mode with m ^ 0. 
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V. NUMERICAL METHOD 
A. The angular part: search of fj, 

The angular part of the perturbation equation, (|2T|) . can be written as 

(1 - {2 + 0^(1 + x)\} Q"{x) + - b - {D - l)x + y (1 + x) {D - 2, - {D + l)x} Q'{x) 

'fjL a'^uj'^ix-l) m2(l + a2A) l{l + D - h) 



(44) 



2 2(2 + a2(l + a;)A) x-l l + x 

in terms of the variable x — cos(20). 

The equation (j44p has four regular singular points: x = ±1, x = ^, and x = oo, where 



Q{x) = 0, 



A = -4- (45) 



i?2 1+^' ij2 

If we introduce the rescaled mode function y{z) by 

Q{X) = (1 - X)l"l/2(1 + ,;)'/2(^ „ ^)±flW2j^(^)^ (46) 

with a; = 2z — 1, then this mode equation can be put into the standard form of the Heun's eauation[35l [3^ . 



2 



70 e 
- + - + 

Z Z — 1 Z — Zq 



d aPz — q 



dz z{z — l)(z — Zq) 



y{z) = 0, a + /3 + l = 7 + 5 + e, (47) 



where 



a^^{l + \m\±Ruj), (3^^{l + \m\+D-l±Ruj), (48a) 

-f^^{2l + D-3), S = l + \m\, e = l±Ruj, (48b) 

9 = + ^ + ^ {{1 + + \m\+D-3)-f,}, (48c) 

^o=^^=^. (48d) 

2 



Note that since 1 + Xa^ must be positive, zq > 1. Due to the symmetry of the equation (|44p with respect to the 
transformation R —^ —R, both signs in the above equations lead to the same eigenvalues. 

The boundary condition on y{z) (0 < z < 1) can be determined as follows. First, if we put the equation ((44l) into 
the eigenvalue equation form iiQ = LQ and require that the operator L be self-adjoint when is real, we find that 
the norm of Q should be defined as 



7V(Q)oc j dxil + xf^-''^'^ 



\Q\\ (49) 



Next, it is easy to see that the fundamental solutions to behave as Q ~ (2; + l)'^^,(a; + l)~('+-^~^)/2 around 
X = —1 and as Q « (1 — a;)*!'"!/^ around a; = 1, respectively. Among these solutions, those that behave as (a; + 1)'/^ 
around a: = — 1 and as (1 — a:)!™!/^ around a: = 1 (m 7^ 0) become normalizable with respect to N{Q). These conditions 
are identical to the condition that y{z) is regular at z = and z = I. 

Now, under this boundary condition, the eigenvalue fi can be determined as follows. First, following the method of 
expansion of the Heun's equation proposed in (30| . 

y{z) = Y.c,u,{z), u,{z) = F{-j,3+p-r, z) - {-iy^^-^P^''"''-'-'\2z - 1), 

we find the three-term recurrence relation for the coefficients Cj : 

Cj+iaj + CjPj + Cj_i7j = 0, 
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with 



(j + +P - a + + P - /? + + S) 



(50) 



= (2j+p+l)(2,+p-l) ^0J{J+P)-,, (51) 

(.7 + a - l)(j + /? - l)(j + 7 - +P - 1) 



where we defined the quantity 

p = j + S-l = a + l3-e = l + \% 

z 

For a given uj, eigenvalues /i can be found through the equation that can be expressed in terms of an infinite 
continued fraction as 

aoji ai72 f . 

/3o-^^.... (53) 

Since the continued fraction is convergent (30| . we are able to find the right-hand side of the equation (j53p with any 
desired precision by limiting the depth of the infinite continued fraction by some large value, always ensuring that an 
increase in this value does not change the final results [3l|. In order to find the j"* root we invert ((53)) j times 

^ aj-i7j aj-2lj-i ao7i ^ aj7j+i «j+i7j+2 ^^^^ 

We shall use this numerical procedure in order to find /i for each given uj. The procedure is convergent and therefore 
accurate for all to and a, and docs not contain any approximation. 

It is also useful to know the non-rotating limit as a start point for numerical search. For that case, a ^ 0, and the 
quantities zq and q are dominant. The equation (|54p is then reduced to Pj = 0, which reads 

q = ~zo3{j+p). j = 0,1,2,3... (55) 

This equation can be easily solved with respect to /i as 

pi={l+ \m\){l + |m| + £> - 3) + A]{j+p) = (2j ^ I ^ \m\){l + 2j 4- |m| + D - 3), (56) 

which does not depend on lo and R in the non-rotation limit. 



B. The radial part: the search of lo 



In terms of P{r) defined by 



P(r) =r-^/2+2(^2^^2)-l/2p(^) 



the radial equation (|22p can be put into the form 



d2p 



2Mam 



(r^ + a^)^r-^ (r^ + a^)^ 



Mr) 



U{r) 



P = 0, 



where 



A(r) 



is the tortoise coordinate, and U{r) is given by 



2M 



V 2 /V 2 / r2 r2 -I- 
{(D-2)r2-h(7:i-4)a2}2_8a2r2 M 



2(r2 a2)2 



r.r>-3 ■ 



(57) 
(58) 



(59) 
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Here, note that the background black hole solution has three parameters M, a and A. From this point, we use the 
horizon radius r^j of the black hole in place of the black hole mass M that is related to by 

2M = r^-\rl+a^){l-\rl). (60) 

We often use the curvature radius R instead of the cosmological constant A = ^, which should be in the range 

a < R due to the regularity condition (see Fig. [1]). 

Now, let us specify the asymptotic conditions on P(r) in order to determine the frequency uj. In the asymptotically 
flat case, these conditions for unstable modes and for quasi-normal modes are formally identical and given by the 
ingoing condition at horizon and the outgoing condition at spatial infinity pretending co to be real. The only difference 
for these two types of modes are just the sign of the imaginary part of the frequency lo. In contrast, for AdS black 
holes, some subtlety arises concerning the boundary condition at the spatial infinity. 

First, at the horizon, we have 



Pir) - (r_^,^)±-*('-h+-^)/A'('-h)^ 



where is the shifted frequency defined in (|29|) . Hence, the quasinormal boundary conditions at the horizon are 
given by 

P{r) = (r - ^^)— 4'-^+a^)/A'(r,)(^^ _^ _ ^^^y (g^) 

Next, at the spatial infinity, the two linear independent solutions of P{r) are 

Pi(r) ^ P2{r) ^ r(^-2)/2^ (g2) 

As far as the instability is considered, it is natural to look for a mode normalizable with respect to the norm that 
makes the equation ([55]) self-adjoint. This requirement is equivalent to the condition that r^^/^P is bounded at 
infinity. Hence, it is natural to require the Dirichlet-type boundary condition at infinity: 

P{r ^ oo) ocr-"/^. (63) 

We will also use this condition when we look for quasi-normal modes. Notice that the choice of this specific fall-off 
is stipulated by the interpretation of QNMs by the AdS/CFT correspondence, and was used previously in [1,0 for 
non-rotating asymptotically AdS black holes. 

Now, let us explain the numerical procedure used to determine lo by solving (j58p under these boundary conditions. 
For convenience, let us introduce the new function 

y{r)=(l-^) Pir). (64) 



If P{r) satisfies the quasi-normal boundary conditions, y(r) is regular at the event horizon. Since the function y(r) 
satisfies the linear equation, we fix its scale as 

y{rh) = 1. (65) 
After the scale is fixed, we expand y{r) near the event horizon as 

yir) = l + ±C,{l-'f)\o{l-'-^y . (66) 

This expansion allows to find y{r) and y'{r) in some region < r < ri with a high precision, if is close enough to 
r/j. Values y(rj) and y'(rj) are used as a boundary condition for the equation (|58p. Since /i(w) is found numerically, 
we solve the equation (I58p for each uj using the ND Solve built-in function in Mathematica for ri < r < rj , where 
rf > rh. 

We have used the following values for the fundamental modes: p = 20, rf = 5000r?i, r^/r^ = 1-1- r^/rf, making 
sure, that the result does not change if p or r/ increases. For higher overtones, larger values of r/ are necessary in 
order to achieve the convergence of the whole procedure. 

In the general case, the behavior of P{r) at infinity is a superposition of the two solutions ([5^ for the Dirichlct 
and non-Dirichlet boundary conditions, Pd{t) and Pn{r): 

P{r) = ZoPoir) + ZMr), (67) 
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where PDij) satisfies the quasi- normal boundary condition (|63p. If uj is the quasi- normal frequency, the corresponding 
solution must satisfy the boundary conditions (j63p at the spatial infinity and, thereby, Z„ = 0. 

Thus, our numerical procedure is the following. We integrate the equation (|58p numerically, using the boundary 
conditions at r = r^, which correspond to the purely in-going wave at the event horizon. At large distance, we compare 
the obtained function '^{rf) with (|67p and find, thereby, the coefficients Zd and Z„ for any given value of w. The 
quasi-normal modes correspond to the roots of the equation 

Zn{Lj) = 0. (68) 

Here, note that as to is generally complex, Z„(aj) is complex as well. 

VI. QUASINORMAL MODE ANALYSIS 

From here and on, we shall imply that 

to = Re{uj) +ilm{uj), (69) 

so that Im (lu) > corresponds to an unstable (growing) mode, and Im (cu) < corresponds to a stable (damped) 
mode. 

A. Stable modes 

Let us note that if the effective potential is not positive definite, it does not necessarily mean instability. Thus for 
scalar-type gravitational perturbations of the D — 5 Schwarzschild-anti-de Sitter black holes, the effective potential 
approaches — oo at r* =0 (see for instance Fig. 1 in [8(]), yet this does not lead to instability as the quasinormal 
modes are damped for that case [1, Q . 

The stability is not proved analytically for R > a > rf/R. Thus, 1 > a > (in the units i? = 1) is the region 
suspected for instability for nonzero m, and we shall search quasinormal modes in this region for different values of 
the quantum numbers: I = 2,3,..., m — 1,2,..., j = 0,1,2,.... Note that the results for the opposite rotation case 
with a —a can be obtained just by reversing the sign of m in those for a > 0. 

We shall also distinguish the following three regimes of asymptotically AdS black holes: 1) large black holes with 
rh ^ R, 2) intermediate black holes with rt = R, and 3) small black holes with rh <C R (see Fig. [Ij. This 
classification, although conventional, is important, because the properties of the quasinormal spectrum are totally 
different for small and large black holes, while an intermediate regime is a kind of transition region between the other 
two. 

For example, for large asymptotically AdS non-rotating black holes, the ratio of the frequency w to the horizon radius 
r/j is roughly independent of rh for the tensor and vector gravitational quasinormal modes [8|. We are convinced by 
numerical computations that the same result holds for rotating black holes. Indeed, for instance, for the fundamental 
mode in the D = 7, a = 0.7 case, we obtained the following values of w/r/^: 5.10783 — 2.62731i, 5.07295 — 2.6235i, 
5.05571 - 2.62124i, 5.04554 - 2.61974i for rh = 10, 14, 18, 22, respectively. This indicates the convergence of the 
ratio Lo/rh to a constant in the limit of large black holes. 

As a — > 0, and the size of the AdS black holes goes zero rh — > 0, the quasinormal modes should approach their pure 
AdS values, i.e., the normal modes of AdS space-time without a black holes, as it was proved in [33|] for the scalar 
field and in ^35| for the gravitational field in four spacetime dimensions. In our higher-dimensional case, we have also 
checked this to confirm that in the limit a — and r?i — > 0, the quasinormal frequencies approach their normal AdS 
values for tensor-type gravitational perturbations 

LOnR 2n + D + 21 + m + j - I, rh ^ 0. (70) 

For non-zero rotations, the quasinormal frequency in the limit rh — > is characterized by that in the pure AdS case 
again, but with a finite correction F depending on a: 

LUnR^2n + D + 2l + m + j-l + F{DJ,m,j,a), r/,. ^ 0. (71) 

We can see representative examples of the dependence of the first few QNMs on a for a fixed rh in Fig. [51 In this 
figure, we also observe that higher overtones are naturally more stable, as they damp quicker. 

One should note that here we have considered the limit of rh keeping R constant (R — 1). This is different 
from another small black hole limit for which i? — > oo keeping rh and a finite. The latter limit corresponds to Kerr 
black holes in the asymptotically flat background. 
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a/R 



FIG. 2: Real (left) and imaginary (right) parts of the first five modes (the fundamental mode and four overtones) for m — 10, 
rh = 0.5, D = 7 

Re(w)R 




FIG. 3: Real (left) and imaginary (right) parts of the fundamental modes, for Vh = 0.5 as a function of a, m = 0, 1, ...16, D — 7. 

Before going to the arguments on instability, let us discuss the dependence of quasinormal modes on the parameters 
D, TO, n, rh and a. Figure [3] shows representative examples of the dependence of the fundamental modes on m. One 
can see that the spacing in the frequency uj between adjacent TO-modes at large to becomes independent of m as 

Re (a>„,+i) - Re M ^ R-^ + O ( - ) , Im (cj„,+i) - Im (cj,„) = O (-] < 0. (72) 

\m / \ TO / 

Thus, a mode with higher to has the negative imaginary part with a larger absolute value, so is more stable. This is 
rather unexpected because the increase of to makes the negative dip of the potential deeper as can be seen from (|59p . 
One possible explanation is that because the eigenvalue fi is approximately given by m? at large to as (j56p suggests, 
the increase of to also increases the height of potential barrier, which supercedes the former effect. 

In order to check the formula ([72| for other values of the background parameters, we have considered a large black 
hole rfi — lOR and calculated ojm for to — 1, 2, 3, ...600. As shown in Fig. [U the high-m spacing behaves in the 
same way as above in the units R = 1. This is qualitatively similar to the high-multipole frequency behavior for 
4-dimensional Schwarzschild black holes, which is known analytically. These results suggest that the formula ([7^ 
holds for all values of 13, and a. 



B. Instability 

The first step in the search for unstable modes is to check that there occurs no instability in the non-superradiant 
range of the spectrum. For this, we have performed extensive calculations of |Z„(a;)| in a wide range of complex 
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FIG. 4: Spacing between real (left) and imaginary (right) parts of high m- modes {D — 7, rh = 10, R = 1). In the right panel, 
Im (curn+i) — Im {ujm) is actually negative and its absolute value is plotted. 




FIG. 5: Absolute values of Z„(a;) for m = 10, rn = 0.5, a = 0.8 (left) and a = -0.8 (right). 



values for lu with a positive imaginary part, which corresponds to instabiUty, and found no zeros of |Z„(cj)| if the 
superradiance condition 

Re (w) < mflh (73) 

is not satisfied (see Fig. [5] for representative examples). 

In order to get an initial hint on where to look for superradiant unstable modes, as the next step, we used the 
following estimation. Since the suspected region of instability is 1 > a > in units R = 1, we are dealing with 



(Re(w)-iiixIlh)xR, In<w)RxlO" (Re{w)-nKnh)xR, Iirtw)Rxlo" 




FIG. 6: Fundamental quasinormal modes m = 2(left) and m = 3(right), ru 
ii(Re(tj) — m£Lh), red line is oc RIm{(jj). 



= 0.1 as a function of a {D — 7). Blue line is the 
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FIG. 7: Fundamental quasinormal modes m — 4(left) and m — 5(right), = 0.1 as a function of a (D — 7). Blue line is the 
R(Re{uj) — mQ,h), red line is oc RIm{Lu). 



(Re(w)-mxnh)xR, Ini(w)RxlO'^ (Re(w)-mxnh )xR, In<w)RxlO 




FIG. 8: Fundamental quasinormal modes m = 2(left) and m = 3(right), D = 8, = 0.1 as a function of a. Blue line is the 
R(Ile{uj) — mQh), red line is oc Rlm(Lo). 



small black holes including those for which the QNMs are close to the pure AdS normal modes, that is, Re (oj) w 
2ji + D + 21 + m + j — 1. For these black holes, the superradiance inequality Re (w) < m^lh leads to the inequality 
2n + D + 21 + m + j — 1 < mflh, which gives us an initial guess in the numerical search for unstable modes with small 
a. In particular, from this inequality and ^Ih = a(l + ''"h)/{^h + o^) > I7 we find that the smaller the black hole is, 
the smaller the minimal value of m for superradiance is. Thus, for instance for D = 7, rn = 0.5, I = 2 and j = 0, the 
superradiance occurs only for m > 32, while for rh = 0.1 and the same D,l and j, the superradiance starts already 
at TO > 2. The above inequality also implies that for higher D and the same fixed rh, the minimal value of m for 
the superradiant instability slowly increases with D. Thus, for example, for D = 8, = 0.1, I = 2 and j = 0, the 
minimal azimuthal number for the instability is to = 3, while to = 2 is not in the superradiant regime (see Fig. [5|) 

(Re(w)-mxnh)xR, InXw)RxlO'-' (Re(w)-mxnh )xR, In<w)Rxl0" 




a 



FIG. 9: Fundamental quasinormal modes m = 4(left) and m = 5(right), D = 8, = 0.1 as a function of a. Blue line is the 
-R(Re(Lj) — mQh), red line is cx RIm{Lo). 
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(Re(w)-mxnh)xR, Irn(w)RxlO 
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(Re(w)-mxnh)xR, Ini(w)RxlO 



-10 



-20 



20 



10 




-40 



a 



FIG. 10: Fundamental quasinormal modes m — 5, D — 10, — O.l(left) and D — 7, m — 2, — 0.02(right) as a function of 
a. Blue line is the R{Ke (uj) — mSlh), red line is cx Rim (lj). 

Now, we summarize the basic features of the results of our numerical calculations. 

1) The Imw is positive (unstable) only for modes in the superradiant regime (see Figs. [&l llOp . 

2) The instability growth rate is tiny and of order 10~^^ or less in units R ~ 1. 

3) The maximal value of the instability growth rate, F, decreases when D increases, at least in the calculated cases 
with D in the range 7 < £> < 10. For example, for D ^7,rh= 0.1 and a w 0.21 in units i? = 1, F = 1.8 x 10"", 
while for higher D, the instability growth is smaller (see Fig. [51 [51 [TU)) . 

4) The unstable range of a increases when we increase m. Black holes becomes unstable for slower rotations with 
the increase of the azimuthal number m. 

The last observation together with the numerical formula (|72p allows us to guess an analytic expression for the 
region of instability. In fact, the superradiance inequality (j73p for large m reads 



Since the instability region in the parameter space expands with the increase of m, the m — > oo limit gives us the 
exact formula for the unstable region 



For any values of a, r^ and R in this region, we can find a value of m that satisfies the superradiance inequality, and 
therefore, the corresponding perturbations are unstable. 

Taking into account the regularity condition a < R and the equation ()43|) . we conclude that the parameter region 
of instability for rotating AdS black holes against tensor-type gravitational perturbations is given by the inequality 



This corresponds to the white (top-left) region on Fig. [TJ 

The tiny instability growth takes place also for a massive scalar field around an asymptotically flat Kerr black 
hole. Thus, it is probable that a small growth rate is generic, although it is a poorly-understood feature of the 
superradiant instability of black holes. Note that such a tiny classical instability of small asymptotically AdS black 
holes will not probably affect the dynamic of small black holes, for which the quantum thermodynamic instability and 
Hawking evaporation are much more effective and violent processes. Moreover, even pure classically, the superradiant 
instability of a black hole in the AdS box does not necessarily mean that the black hole disappears or collapses to 
singularity. At least near the border of the unstable parameter region (the white top- left region on Fig. [T]), it would 
be natural for a black hole to diminish its angular momentum until it reaches a stable state, which is still described 
by a Myers-Perry- AdS solution with smaller a. The complete and definite answer on the evolution and destiny of 
rotating black holes in AdS space should be given by full non-linear analysis of perturbations. 




nhR>i. 



(74) 



rl<aR< 



(75) 
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VII. CONCLUSION 



In the present paper, we have derived the master equation for the tensor- type gravitational perturbations of a 
simply rotating Myers-Perry black hole with the A-term in the spacetime dimension D > 7. We have analytically 
proved the stability of black holes obeying inequality a < rf^/R, while for small rotating AdS black holes violating 
this condition, we have found unstable modes. We further found that the parameter range for instability exactly 
coincides with that for superradiance instability. This result indicates that one should be careful when one utilizes 
a stationarity of the black hole metric as an implicit background for propagation of fields, for example when one 
calculates the grey-body factors (which implicitly provide the emission probabilities for the Hawking radiation) or 
quasinormal modes of higher dimensional black holes [13] . 

We have also found a strong numerical evidence that the parameter region for unstable black holes is identical to 
the region Rflh > 1 or equivalently a > r\/ R in which the analytical stability proof does not hold. In addition, we 
have analyzed the spectrum of the damped quasinormal modes in the stable sector and found some general tendency 
in the dependence of the complex frequency on the mode parameters and the black hole parameters. 

Although we have found some parameter region of instability for rotating AdS black holes here, it is possible that 
the complete parameter region of instability is larger than the one found here, because we have considered only tensor- 
type gravitational perturbations. In particular, our analysis is powerless in studying the stability of rotating black 
holes in five or six dimensions. In order to investigate the stability of these black holes, it is clear that the scalar-type 
and vector-type gravitational perturbations have to be studied. Unfortunately, these are at present out of our reach 
because perturbation equations in those cases have not been reduced to a set of tractable ODEs by the separation of 
variables except for some exceptional cohomogeneity one cases such as the special five-dimensional Myers- Perry black 
holes with ai = 02 [stI- 

In spite of this restriction, there still exist a lot of interesting problems that can be handled by the present formu- 
lation. Some examples are the investigations of the tensor-type gravitational quasinormal modes for simply rotating 
black holes in the asymptotically flat case and in the asymptotically de Sitter case. The above formulation can be 
also used to calculate the classical scattering amplitudes as well as the probability of emissions (due-to the Hawking 
radiation) of tensor-type gravitons from a simply rotating black hole. 
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